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The unphysical solutions of the periodic Anderson model 
obtained by H. Keiter and T. Lenders [Europhys. Lett., 49, 
801 (2000)] in dynamical mean-field theory (DMFT) are 
shown to result from the author's restricted choice of the func- 
tional form of the solution, leading to a violation of the an- 
alytic properties of the exact solution. By contrast, iterative 
solutions of the self-consistency condition within the DMFT 
obtained by techniques which preserve the correct analytic 
properties of the exact solution (e.g., quantum Monte-Carlo 
simulations or the numerical renormalization group) always 
lead to physical solutions. 



In a recent paper, Keiter and Lenders [Q discussed the 
structure and the solutions of the self-consistency condi- 
tion of the dynamical mean-field theory (DMFT) for the 
periodic Anderson model at t7 = oo. They raised two in- 
dependent points: (i), they showed analytically that on 
certain assumptions for the structure of the solution the 
self-consistency condition can have multiple solutions, 
which merge in the upper half of the complex plane, and 
(ii), they argued that calculations of 1/d corrections (d: 
spatial dimension of the lattice) to the self-consistency 
condition of the DMFT are problematic, since the limit 
d ^ oo does not commute with the limit n ^ oo of the 
length n of loops embedded in the lattice. 

In the following both points will be addressed. The 
main part of the paper will be concerned with a detailed 
analysis of the first issue. In particular, we will show 
that the conclusions of Keiter and Lenders |^] relating to 
this point are unwarranted, i.e., are due to the deficiency 
of the method used by these authors to solve (or simu- 
late the solution of) the single-impurity Anderson model 
(SIAM). 



A. Structure of the solution of the self-consistency 
condition 

Multiple solutions of a nonlinear set of equations are 
not unusual. In particular, most multiple solutions of 
nonlinear equations for propagators are unphysical, i.e., 
have poles in the upper half-plane or lead to an incor- 
rect asymptotic behavior (or both). In this respect the 
findings by Keiter and Lenders |0] that under certain con- 
ditions there exists a branch point in the upper half-plane 
where an unphysical and the physical solution to the self- 
consistency condition merge, is disturbing. If true in gen- 



eral it would raise doubts about the uniqueness of numer- 
ically obtained iterative solutions to the self-consistency 
condition, and would make the physical relevance of re- 
sults obtained within the DMFT questionable. 

Here we will show that, in fact, the self-consistency 
condition for the periodic Anderson model (PAM) |^ al- 
ways leads to a unique physical result unless the func- 
tional space for the solutions is artificially restricted as 
done in ref . [0 . In standard notation the PAM reads 



PAM 



kcr 



(1) 



Within the DMFT, the nontrivial part of the partition 
function is represented by an effective local action j^] 

= - E / drdr' fl{T)g-\T - r')f,{r') (2) 
+[/ / dr/|(r)/^(T)/t(r)/^(r) 



where 



Q[z)=[z-es-V^g{z 



(3) 



and the auxiliary propagator Q{z) describes an artificial 
system where the /-states at the site under consideration 
have been removed. Note that the action (|^) defines an 
effective SIAM with band-electron Green function g{z) 
which has to be determined self-consistently. For the 
PAM, (|l|), Keiter and Lenders derived the following ex- 
pression ||l[: 



-E- 

N ^ 1 



-g^{z)T{z) + T{z)gizy 



(4) 



Here gk{z) — l/(z — ek) is the propagator of the non- 
interacting conduction band system and T[z) — V'^F{z), 
with F{z) as the local /-electron Green function. The 
latter relation implies that T{z) cannot be chosen at lib- 
erty. Eq. (^ is not sufficient to solve the problem, since 
F{z) is unknown. The quantity F{z) has to be deter- 
mined from an effective SIAM, where the propagator of 
the host electrons is given by g{z). It should be noted 
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that the relation is not only valid ior U = oo, but 
for all values of U, including U = 0. Moreover, it is 
important to mention that in Ref. the artificial host 
propagator Giz) appearing in (|) is denoted by G{z), 
while conventionally, and also in this paper, G{z) refers 
to the propagator of the physical band states. 

Let us briefly comment on the validity of Equ. (Q) 
for U = 0. At first glance it seems rather strange that 
the solution to a non interacting system should be given 
by a complicated, non-linear self-consistency condition. 
The reason, however, is quite simple. First, Equ. (Q) 
is set up in real space, while the solution to (Q) for 
[/ = is naturally obtained in reciprocal space. Second, 
as mentioned before, the quantity G{z) is not a physi- 
cal propagator but an artificial auxiliary quantity. It is 
perfectly valid, although not necessarily reasonable, to 
rewrite the solution of the PAM (]l|) in terms of this aux- 
iliary local propagator even for U — 0, thus obtaining 
the complicated expression (jj). For the physical prop- 
agator G{z), on the other hand, one eventually finds. 



inserting T{z) = V'^ 



6/ — y^tJ(z)] , a simple and 



sensible solution (cf. Equ. (|7|) below). It must be noted 
that for U = this procedure is exact independent of 
dimension provided the hybridization in (|^) is local. The 
DMFT or limit d = oo enters at a different stage, namely 
through the locality of T(z) even in the presence of the 
two-particle interaction U. 

To examine the analytic properties of (|^) under the 
condition that T{z) is given by the solution to an effec- 
tive SIAM some simple algebraic transformations have 
to be performed. Let us assume that T(z) and G{z) are 
solutions to (0). Then we may rewrite (0) as 
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l + T{z)G{z)N^' 



i+T{z)g(z) 



(5) 



For the SIAM, the host-electron propagator of the inter- 
acting system is obtained from the host propagator of the 
system without /-states by 



Giz) = Giz)+Giz)Tiz)Giz) 



(6) 



We can rewrite (^) in terms of the physical host-electron 
propagator as 



G{z) 



1 



1 



i+T(z)g{z) 
and thereby identify the self-energy 

T(z) 



l + T{z)G{z) 



(7) 



(8) 



of the host electrons. Note that for U — the self en- 
ergy (^) reduces to the well-known expression 'E^'^^z) — 
V"^ l\z — tj\. Since the k-dependence in is only due to 
ek, the k-sum can be converted to an energy integral 



-E- 



(9) 



with a positive-definite density of states (DOS) p'-°^(e). 
The relation (|^) is simply a Hilbert transform; therefore 
one can make use of its analytical properties to show ||^ 



G{z) 



(10) 



SM(z) -r(z- S]M(z))' 

with a Herglotz-analytic function r(C), i.e. 

sign(3mr(C)) — — sign(3mC) 

and r(C) ^ 1/C for |C| ^ co- Furthermore, r(C) is, in 
principle, a known function which is determined by the 
structure of the DOS alone. Finally, G^z) can be obtained 
from the definition of G{z)^^ — Giz)^^ ~ I](^^(z) as 



G{z)-' = G{z)-' + E(^)(z) = z-r{z- 



T{z) 



\ + T{z)G{z)' ■ 
(11) 



It should be stressed that (|llD is not merely a nonlinear 
equation determining G{z) but, rather, a functional equa- 
tion, since T{z) is a functional of G{z). The functional 
dependence is determined by the solution of a nontriv- 
ial problem, namely the SIAM. It is thus not sufficient 
merel y to study the possible roots of the nonlinear equa- 
tion ([ll|), but one has to discuss the structure of this 
equation under the constraint that T{z) must obey the 
basic properties of the solution to the SIAM (||). Since 
there does not exist an analytic solution to the SIAM in 
closed form this constraint can be incorporated currently 
only within an iterative scheme, even if T{C,) is known in 
closed form as, for example, for the Bethe lattice or the 
toy DOS used in [|. 

Let us therefore assume that we knew the exact solu- 
tion of the SIAM, and study the consequences of an it- 
erative procedure to solve the self-consistency condition. 
Starting in the n-th step of this iterative procedure with 
a physical C/'-"-'(z) with sign (3mCj(")(z)) = — sign(3m2;) 
and G^''\z) - 1/z (|z| oo), the corresponding propa- 
gator of the /-states, F("^(z), and hence T("^(z), will 
have the same properties, i.e., the right-hand side of 
(O) will definitely have the correct asymptotics, because 
tF0(^)5(»)(2) _ 1/^2 ^Yms 



1 -f T(")(z)e;(")(z) 



1/z 



as l^l — » oo. More importantly, however, the exact solu- 
tion to the effective SIAM (||) has to satisfy the condition 



sign 



(T(")(z)-i-f ^(")(z))] -+sign(3TOz) (12) 



which is necessary to ensure that I](^^(z), (^), has the 
correct analytic properties, i.e., sign (SmS*^"^^ (z)) = 
— sign (3mz). Under those conditions, the iterative pro- 
cedure always stays within the subspace of physical so- 
lutions. Here it must be stressed once more that one is 
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not free to choose T'^'^^z) = V^F'-'^^z), or even choose a 
particular functional form for r^"^(z), since this form is 
entirely and uniquely determined by the solution of the 
effective SIAM with a fixed host Green function Q^"''{z). 
Especially condition (12) turns out to be very susceptible 
to approximations and is easily violated. Such a violation 
would manifest itself in an indefinite sign of 3mE(z) and 
would thus lead to the crossing of a branch point into 
the upper half plane. Note also that in the course of 
the iteration, the quanitities entering on the right-hand 
side of ( pi] ) are those of the n-th step, producing the host 
Green function CJ^""'"^^(z) for the n + l-th step and, unless 
self-consistency has been reached, ^/^""'"^^(z) ^ C/'-"-'(z) in 
general 0. 

Hence we find that the branch point in the upper half- 
plane observed by Keiter and Lenders Q is merely due 
to the author's restricted choice of the functional form 
of T{z) which leads to a violation of (p^. Namely, mo- 
tivated by results for the standard SIAM these authors 
chose for T{z) the following ansatz j|] 



self-consistency condition and not of the self-consistency 
condition itself. 



T{z) = 



Ok 



iTi. 



(13) 



This form represents the strong-coupling regime of the 
SIAM, i.e. -e/, Cf+U > 7rp(o)(0)F2. The original single- 
particle level at e/ is broadened by the coupling to the 
conduction bath. In addition, the sharp Abrikosov-Suhl 
resonance which characterizes the virtual bound state 
generated by spin-flip scattering appears at the Fermi 
energy. A typical result will have the form shown in 
Fig. I^a. In Fig. |] we assumed a Gaussian free DOS 
/o(°He) = exp(-e2)/0F and set = 0.2, £/ = -1. Let 
us compare this result with a fully self-consistent DMFT 
calculation for the PAM with the same parameters 
The result is shown in Fig. |l]b (full line). It is obvious 
from the inset to Fig. |^b that the true solution of the 
PAM does not have the simple structure given by ( p^ ) 
assumed in ref. ||l[. In particular, one does not observe a 
single peak at the Fermi energy, but two peaks separated 
by a gap. 

Consequently, the conclusion following equ. (10) in ref. 
is unwarranted: the observed breakdown, i.e. the ap- 
pearance of a branch point in the upper half plane, does 
not result from the fact that in this regime physically in- 
terpretable single-valued solutions cannot be found, but 
merely from the deficiency of the method used to solve 
(or simulate the solutions of) the SIAM. In particular, 
one of the authors jl^ and F. Anders jlll observed that 
approximations based on the resolvent perturbation the- 
ory ]l^ , like the non-crossing approximation (NCA) or 
improved versions of it, always lead to spectra similar to 
( p^ ) and are not able to reproduce more complicated fea- 
tures like the double-peak structure at the Fermi energy 
appearing in Fig. [^b. It is therefore not surprising that 
the application of these techniques in solving the DMFT 
for the PAM leads to the breakdown described in ref. : 
it is a problem generated by the method used to solve the 
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FIG. 1. Imaginary part of T{Lu + iS) for the PAM. (a) result 
according to (n^ as used in ref. (b) result from a full 
DMFT calculation with the NRG g. In the inset in (b) 
the two curves are compared for energies close to the Fermi 
energy. 

Our arguments do not imply that the self-consistency 
condition (|l^) must necessarily have only one unique so- 
lution for a given set of model parameters. Such a con- 
clusion would require a more detailed knowledge of the 
solutions of the SIAM (especially, whether there exists a 
one-to-one correspondence between Giz) and F{z)), and 
of the extent of the basin of attraction for a particular 
solution. Since there does not exist an exact analyti- 
cal solution to the SIAM for an arbitrary host Green 
function G{z), statements about these properties of the 
self-consistency condition are out of reach. Indeed for 
the Hubbard model such a situation (namely, two physi- 
cally valid, coexisting solutions) occurs close to the Mott- 
Hubbard metal- insulator transition Hi . 



B. 1/d corrections to the DMFT 

In DMFT a fermionic lattice model is mapped onto a 
SIAM, complemented by a self-consistency condition 
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This mapping is rigorous in d = oo, a fact which is not 
disputed in ref. 0|. 

The possibihties and Umitations of systematic \/d- 
expansions around d = oo have aheady been explored by 
a number of authors ||l3|-[l9t. Whether l/c?-corrections 
to DMFT are able to improve the theory depends on 
the quantity under investigation, on the type of 1/d- 
expansion, and also on the physics one would like to de- 
scribe. For example, while 1/d-corrections include non- 
local short-range correlations they are unable to cap- 
ture certain long-range correlations. Furthermore, 1/d- 
corrections, when truncated at finite order, may violate 
basic analytical properties of dynamical quantities [ p^ . 
All this is not unexpected. Indeed, analogous problems 
arise in 1 /c?-expansions around the Weiss mean- field the- 
ory for the Heisenberg model and the coherent poten- 
tial approximation for the Anderson disorder model . 
We do not see why the peculiarities of large loops in 
the large d limit emphasized by Keiter and Lenders 0] 
should make expansions around DMFT particularly dif- 
ficult or ambiguous. Since even the Green function of 
the non-interacting system can be expanded in loops, it 
is clear that these peculiarities are in one-to-one corre- 
spondence with the special behavior of the bare DOS in 
the limit d — *■ oo (e.g., the absence of van Hove singular- 
ities) |2^,|2^ . Potential problems with the 1 /c?-expansion 
of the bare DOS can be avoided by inserting the actual 
DOS of the finite-dimensional system into the DMFT 
equations. Then only corrections to the local approxima- 
tion of irreducible vertices underlying the DMFT need 
to be considered, as has been done already by several 
authors |||,|8|j22||2|l . 



In summary, we showed that the iterative solution 
of the self-consistency condition within the DMFT al- 
ways lead to solutions that stay in the physical subspace. 
Therefore, the unphysical branch cuts observed by Keiter 
and Lenders ||] are due to the method used by these 
authors to solve the SIAM which violates fundamental 
analytic properties of the exact solution. Similar vio- 
lations are caused by the NCA and other approxima- 
tions based on the resolvent perturbation theory . In 
fact, this also appears to be the reason why it was not 
possible up to now to obtain meaningful results for the 
PAM using NCA and related methods, while techniques 
which preserve the correct analytic properties of the ex- 
act solution (e.g., iterated perturbation theory, quantum 
Monte-Carlo simulations or the numerical renormaliza- 
tion group) do not encounter similar difficulties. 

We thank F. Anders, V. Janis, M. Jarrell, H. Keiter 
and D. Otto for useful communication. 
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